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Abstract 

Using data collected at the ^(3770) resonance with the CLEO-c detector at the Cornell e + e" 
storage ring, we present searches for 25 charmless decay modes of the Y>(3770), mostly multibody 
final states. No evidence for charmless decays is found. 
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The -0(3770) is the lowest-mass charmonium resonance above DD threshold jl]. It may 
be the 1 3 D\ state or a mixture of 1 3 D\ and 2 3 Si. Charmless decays of the -0(3770) can shed 
light on S — D mixing, missing ip{2S) decays such as 0(2,5) — > pn, the discrepancy between 
the total and DD cross section at the -0(3770), and rescattering effects contributing to an 
enhanced b — > s penguin amplitude in B meson decays 0. 

The total cross section at the 0(3770) was estimated from older measurements to be 
(7.9 ± 0.6) nb in |2j, which is larger by about 20% than the measured DD cross section 
(6.39 ± O.lOlo^g) n b 0. While the significance of the discrepancy between the total cross 
section and the DD cross section is not large, identifying non-DD decays of the ^(3770) will 
place the discrepancy on a solid footing and shed light on the nature of the ^(3770). The BES 
Collaboration observed ^(3770) — > 7t + tt~ J '/ip with a branching ratio of (0.34±0.14±0.09)% 
0, while CLEO measured a branching ratio of (0.189 ± 0.022±^)% @. This non-DD 
channel contributes approximately 100 keV to the -0(3770) decay width, which motivates 
the search for other exclusive final states. 

Charmless decays of the -0(3770) may provide an avenue to study rescattering effect 
relevant to B meson decays. For example, the rj' exhibits enhanced production in charmless 
inclusive and exclusive B meson decays which is not well understood. If the ■0(3770) decays 
to DD pairs which subsequently re-annihilate into non-charmed final states, a similar effect 
could contribute to enhanced b — > s penguin amplitudes in B meson decays, including modes 
containing an r]', i.e. b — > ccs — > qqs, where q = (u, d, s). 

Although the -0(3770) is believed to be primarily the l 3 Di state of the cc system, its large 
leptonic width indicates mixing with S'-wave states, particularly the ip(2S). By enhancing 
the rate of non-DD decays at the 0(3770), mixing between the 0(3770) and the 0(2S') pro- 
vides an explanation for the anomalously small 0(25') branching fractions to some hadronic 
2-body final states such as pn 0. 

In this Letter, we report results of searches for a wide variety of non-DD modes, including 
final states with and without strangeness and with and without baryons. The modes r/(')37T 
are included since recent predictions exist 0. Modes with baryons are included since any 
observation would be unambiguously a non-DD decay because D mesons are not sufficiently 
massive to decay to baryon pairs. 

The data sample used in this analysis is obtained at the ^(3770) and the nearby contin- 
uum in e + e~ collisions produced by the Cornell Electron Storage Ring (CESR) and acquired 
with the CLEO-c detector. The CLEO-c jf| detector is a modification of the CLEO III de- 
tector 0, in which the silicon-strip vertex detector was replaced by a six-layer all-stereo 
drift chamber. The solid angle coverage for charged and neutral particles is 93% of An. 
The charged particle tracking system, operating in a 1.0 T magnetic field directed along 
the beam axis, achieves a momentum resolution of ~0.6% at p = 1 GeV. The calorimeter 
attains a photon energy resolution of 2.2% at = 1 GeV and 5% at 100 MeV. Two particle 
identification systems, one based on energy loss (dE/dx) in the drift chamber and the other 
a ring- imaging Cherenkov (RICH) detector, are used together to identify 7r ± , K ± and p(p). 
The combined dE / dx-RICR particle identification procedure has efficiencies exceeding 90% 
and misidentification rates below 5% for these three particle species. 

The integrated luminosity (£) of the datasets was measured using e + e~, 77, and p + p~ 
final states [8j]. Event counts were normalized with a Monte Carlo (MC) simulation based on 
the Babayaga event generator combined with GEANT-based detector modeling. The 
data consist of £=55.8 pb _1 on the peak of the -0(3770) and 20.70 pb _1 at the continuum 
(•^/i=3.67 GeV), which is used for background subtraction. The nominal scale factor used 
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to normalize continuum yields to ^(3770) data is / c ° om = 2.55 ± 0.26, and is determined 
from the integrated luminosities of the data sets corrected for an assumed 1/s dependence 
of the cross section where the error is from the uncertainties in the relative luminosity and 
the s dependence of the cross section. The scale factor differs by a factor of 5.2% for each 
power of 1/s. The value of f co used for each mode also corrects for the small difference in 
efficiency between the t/>(3770) and continuum data. 

The analysis strategy and selection criteria are the same as in the CLEO-c analysis 
of exclusive hadronic decays at the ip(2S) Standard requirements are used to select 

charged particles reconstructed in the tracking system and photon candidates in the Csl 
calorimeter. We require tracks of charged particles to have momenta p > 100 MeV and to 
satisfy |cos#| < 0.90, where 9 is the polar angle with respect to the e + direction. Each 
photon candidate satisfies E 1 > 30 MeV and is more than 8 cm away from the projections 
of charged tracks into the calorimeter. Particle identification is used for all charged particle 
candidates. Pions, kaons, and protons must be positively and uniquely identified. That 
is: pion candidates must not satisfy kaon or proton selection criteria, and kaon and proton 
candidates obey similar requirements. Charged particles must not be identified as electrons 
using criteria based on momentum, calorimeter energy deposition, and dE/dx. 

The invariant mass of the decay products from the following particles must lie within lim- 
its determined from MC studies: tt° (120 < M 77 < 150 MeV), 77 (500 < M 77 < 580 MeV, 
530 < M n+n - w o < 565 MeV), u (740 < M n+n - w0 < 820 MeV [760 < M w+w - n o < 800 MeV 
for the upp final state]), (1.00 < M K+K - < 1.04 GeV), and A (1.112 < M pn - < 
1.120 GeV). For tt° — > 77 and 77 — > 77 candidates in events with more than two photons, 
the combination giving a mass closest to the known 7r° or 77 mass is chosen, and a kinemat- 
ically constrained fit to the known parent mass is made 12]. To suppress electromagnetic 
energy deposits in the calorimeter mimicking a 7r° or rj, each electromagnetic shower profile 
is required to be consistent with that of a photon. For 77 — ► 7r + 7r~7r° and uj — > 7r + 7r~7r , the 
7T° is selected as described above, and then combined with all possible combinations of two 
oppositely charged pions choosing the combination that is closest to the known rj(u) mass. 
For A — > p7t~ , a fit of the p and 7r~ trajectories to a common vertex separated from the 
e + e~ interaction ellipsoid is made. Contamination from K s decays is eliminated by particle 
identification and energy conservation requirements. 

Reconstructed events must conserve momentum and energy. The hadrons comprising 
these events each have momentum p$ and combined measured energy E vis . We require the 
measured scaled energy E v[s /E cm be consistent with unity within experimental resolution, 
which varies by final state. We require |£pi|/i? cm < 0.02. Together these requirements 
suppress backgrounds with missing energy or incorrect mass assignments. The experimental 
resolutions are smaller than 1% in scaled energy and 2% in scaled momentum. 

For the final states with four charged tracks and a ir°, an additional cut is applied to 
remove a background of radiative events. When the highest energy photon in an event is 
combined with a low-energy photon candidate, it can imitate a 7r°. We require (^tracks + 
E^)/E cm < 0.995, where E^ is the energy of the highest energy photon. For the final 
states 2(7r + 7r~) and 2(7r + 7r~)7r°, there is a background from (7)7r + 7r~ J/ip arising mostly from 
radiative returns to the ip(2S). This background is vetoed if the recoil mass against the two 
slowest oppositely charged tracks (assumed to be pions), m| low , satisfies 3.15 < m| low < 
3.22 GeV, and/or the invariant mass of the two fastest oppositely charged tracks (assumed 
to be muons unless the dE/dx measurement is consistent with electrons), m| ast , satisfies 
3.05 < m!f * < 3.15 GeV. For the final states 2(7r+7r-)7r°, K + K~7i + 7i~ and K+K-n+n-ji , 
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in order to remove DD background, we exclude events in which the invariant mass of the 
following combinations of particles is consistent with a D°/D ± meson: 7r + 7r~7r°, K ir^, 
71+71", K + K~, K+Ti+71 , K+K-Ti or 7r ± vr°. 

For every final state, a signal selection range in E vis /E cm is determined by Monte Carlo 
simulation, and a sideband selection range is defined to measure background. The signal 
range in E vis / E cm varies between 0.98 — 1.02 and 0.99 — 1.01 depending on the final state. 
Final states with an intermediate rj, u, or must satisfy a scaled energy signal selection 
range requirement identical to the corresponding mode without the intermediate particle, 
and the event yield is determined from signal and sideband selection ranges of the particle 
mass. For example, the scaled energy signal selection range is the same for <pK + K~ and 
K + K~K + K~ . Most modes studied in this Letter have resonant sub- modes, however, only 
narrow resonances are included in this analysis. 

In Fig. H the scaled energy and invariant mass distributions are shown for two typical 
modes: 2(7r + 7r _ )7r° and K + K~7T + 7T~ . Evidence for production of the uj and resonances, 
respectively, is observed in the corresponding mass spectra. The background from DD is 
tiny, and thus is almost invisible. 
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FIG. 1: Distributions for the modes 2(vr + 7r-)vr [(a) and (b)], K+K-ix+tx- [(c) and (d)]. The 
pairs of arrows indicate the signal regions, (a) and (c) The scaled total energy, (b) The 7r + 7r~7r° 
invariant mass in 2(7r + 7r~)-7r . (d) The K + K~ invariant mass in K + K~tt + tt~ . Filled circle with 
error bar: 0(3770) data, solid line: ^(3770) — > 2(7r + 7r _ )7r° or K + K~7r + 7r~ Monte Carlo, dashed 
line: V(3770) -> DD -> 2(7r + 7r")7r or K+K'tt+iv' Monte Carlo, shaded histogram: e+e" -> 7* -> 
2(7r + 7T - )7r° or K + K^-k+tt^ from continuum data. 



Event totals are shown for both the ^(3770) and the continuum in Table HI where 5^(3770) 
(Sco) is the number of events in the signal region and -8^,(3770) (-B co ) the number of events in 
the sideband region in ■0(3770) (continuum) data. Also shown are yields for a DD Monte 
Carlo sample corresponding to 10 times the integrated luminosity of the data: S DI j (B D jj) 
in the signal (sideband) region. Under the assumption that interference between -0(3770) 
decay and continuum production of the same final state is absent, the number of events 
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attributable to each -0(3770) decay mode, N$, is 



Ns — 5^(3770) — -8^,(3770) — fco{S co — B co ) — foD(S DD — B DD ), (1) 

where f co is mode dependent and listed in Table HJ and /ds=0.1 is the scale factor for DD 
MC. Since no statistically significant excess is observed, we obtain a 90% C.L. upper limit 
on the number of events by adding 1.64 times the statistical uncertainty determined from 
the yields on the continuum and at the "0(3770). 

The efficiency, e, for each final state is obtained using a MC simulation [lOj of the CLEO-c 
detector based on the EvtGen event generator ^3|. No initial state radiation is included 
in the Monte Carlo, but final state radiation is accounted for. The efficiencies in Table H] 
include the branching ratios for intermediate final states. 

We correct the number of events by the efficiency e, and normalize to the integrated 
luminosity £ to obtain the cross section: 

N s 

° = -c 

and normalize to the total number of the 0(3770) decays to obtain the branching ratio: 

N s 



B 



e -^(3770) 



The number of the -0(3770) decays is determined by N^ 770 ) = cr tQt -£, where a tot = (7.9±0.6) 
nb is the world average total cross section of the -0(3770) from [2J. 

The systematic uncertainties on the ratio of branching fractions share common contri- 
butions from the integrated luminosity (1.0%), uncertainty in f co (10.0%), trigger efficiency 
(1.0%), and Monte Carlo statistics (2.0%). Other sources vary by channel. We include the 
following contributions for detector performance modeling quality: charged particle track- 
ing (0.7% per track), 7i /r](^ 77) finding (4.4%), A finding (3.0%), ir/K/p identification 
(0.3%/1.3%/1.3% per identified ir/K/p), and scaled energy and mass resolutions (2.0%). 
The systematic uncertainty associated with the sideband background is obtained by coher- 
ently increasing the backgrounds both at the -0(3770) and on the continuum by one statistical 
sigma. Since the background in many modes is small, the Poisson probability for the ob- 
served number of background events to fluctuate up to the 68% C.L. value is calculated 
and interpreted as the uncertainty in the level of background. Many of the modes studied 
have resonant submodes, however the efficiencies do not differ by much. We generate MC 
data with a phase space model, and take a 10.0% uncertainty for decay model dependence. 
In the computation of the branching fractions, a common uncertainty of 7.6% enters due 
to the number of -0(3770) decays arising from the uncertainty in the total cross section of 
the V(3770) 0. We give the significance for each mode and upper limits (including the 
systematic error) for the cross section and branching ratio in Table |U 

In summary we have searched for 25 exclusive multibody hadronic decay modes at the 
-0(3770). No significant signal is observed in any mode. For each mode we give the signif- 
icance, and the upper limit on the cross section and branching fraction at 90% C.L. This 



study together with the -0(25) multibody decay analysis [ll( provide useful information 
about S — D mixing. The cross section deficit remains a puzzle. However the uncertainty in 
the total cross section of the 0>(377O) is large. A fine energy scan over the -0(3770) resonance 
to measure the total cross section would be very valuable. 
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TABLE I: For each final state h the following quantities are given: the number of events in the 
signal region, SW3770), and background from sidebands, -B,/>(3770)> in ^(3770) data; the number of 
events in the signal region, S co , and background from sidebands, B co , in continuum data; the scale 
factor, / co ; the number of events in the signal region, S d q, and background from sidebands, B D f), 
in a DD MC sample corresponding to 10 times the integrated luminosity of the ^(3770) data 
sample; the number of events attributable to -^(3770) decay, Ns, computed according to Eq. ^ the 
significance, in units of standard deviations; the efficiency, e; the cross section upper limit including 
the systematic error (90% C.L.), and the branching ratio upper limit including the systematic error 
(90% C.L.). For 77371", the two decays modes r] — > 77 a and r\ — > 37r b are combined on line r/37r. 
(The sign of the significance indicates an excess/deficit of events). 
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B U.L. 


h 




B co 








S-ip(3770) 










(Pb) 


(10- 4 ) 


2(vr+7r-) 


1471 


28 


2.49 


1 


13 


3411 


90 


-266.5 


-2.5 
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0.1 
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0.0248 


24.0 
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1 





5 
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7.0 
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33 


13 
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25 


3.3 


0.2 


0.1629 


3.2 


4.1 


<t>h 


12 
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32 


15 
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0.0 


0.0863 


3.5 
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634 


18 


1.73 


30 


21 
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0.5 
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23.6 
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3 





2.50 








3 





-4.5 


-0.7 


0.0389 


3.2 


4.1 


ujK+K- 


62 


12 


2.31 
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14 


-15.3 
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2(K + K') 
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11 
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21.7 
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3 
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